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We calculate the current-pressure characteristics of a ballistic pinhole aperture between two vol- 
umes of B-phase superfluid ^He. The most important mechanism contributing to dissipative currents 
in weak links of this type is the process of multiple Andreev reflections. At low biases this process is 
significantly affected by relaxation due to inelastic quasiparticle-quasiparticle collisions. In the nu- 
merical calculations, suppression of the superfluid order parameter at surfaces is taken into account 
self-consistently. When this effect is neglected, the theory may be developed analytically like in the 
case of s-wave superconductors. A comparison with experimental results is presented. 
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I. INTRODUCTION 

Liquid ^He is a strongly interacting system of fermionic 
atoms with nuclear spin 1/2. Its superfluid state be- 
low the critical temperature Tc ~ 1 mK is characterized 
by the creation of a condensate where the atoms form 
Cooper pairs 1]. This is similar to what happens for 
electrons in superconducting metals and, although '^He 
atoms are electrically neutral, many analogues exist be- 
tween the transport properties of the two physical sys- 
tems. For example, in both systems so-called Andreev 
reflection can exist, where quasiparticles are converted 
between particle-like and hole-like branches of the exci- 
tation spectrum by the pairing potential However, 
instead of the singlet s-wave state of conventional super- 
conductors, the pairing state in superfluid ^He exhibits 
spin-triplet p-wave symmetry. This means that the con- 
densate has internal degrees of freedom, resulting in a 
complicated structure for the order parameter, and in 
the existence of multiple superfluid phases. There is also 
no crystal potential to impose symmetry restrictions, as 
in the case of unconventionally paired (d-wave) super- 
conductors. As a result, many of the analogous phe- 
nomena occur in a more complicated form in ^He than 
anywhere else. In this paper we study the properties of 
pressure-biased weak links in superfluid "^He. The weak 
links consist of small apertures in a wall between two 
volumes of the liquid 3], and, as such, are analogous to 
ballistic point contacts between superconducting metals. 
The theory of superconductor point contacts is well de- 
veloped, and thus most of the basic ideas may simply be 
inferred from existing results 0, S S S B 13 ■ 

Most importantly, due to Andreev reflection, there 
are bound quasiparticle states localized at the weak 
link, whose energies are below the bulk gap A [Tol |. 
These sub-gap states are responsible for carrying the 
phase-dependent supercurrents, i.e., the Josephson ef- 
fects 13, ^3- When the contact is biased by a chemi- 
cal potential difference U , the supercurrents oscillate at 
the Josephson frequency 2U/h. Under such a bias, also 
dissipative dc currents will be generated. The most ob- 
vious source of such currents is due to thermally excited 



quasiparticles, but the resulting current is very small at 
low temperatures. However, in the case of a point con- 
tact, the Cooper pairs themselves may participate in the 
flow of a dissipative current. This is because in trans- 
mitting a pair between the two condensates, energy can 
be conserved by transferring the excess energy 2U to the 
bound-state quasiparticles. As a result, large dc currents 
can flow with arbitrarily small biases U <^ A also at low 
temperatures. The underlying process by which this is 
accomplished is known as multiple Andreev reflections 
(MAR) ^9]. In this process, the bound quasiparticles 
are Andreev-reflected several times from the surround- 
ing pairing potentials under the influence of the bias U. 
After two successive retroreflections, a quasiparticle has 
gained the energy 2U, and this corresponds to the dis- 
sipative transmission of one Cooper pair. This coherent 
process is repeated until the bound quasiparticles escape 
to energies above A, or are relaxed by inelastic scat- 
tering. The maximum number n of sub-gap reflections 
is given by nt7 = 2 A. In superconductor weak links 
with non-perfect transparency, MAR can give rise to a 
highly pronounced "subharmonic gap structure" (SGS), 
where the differential conductance is peaked at the biases 
U = 2A/n, with n = 1, 2, 3, ... 0. On the other hand, 
in the limit of very low transparency, a tunnel junction 
is formed, where the sub-gap states and hence MAR are 
completely suppressed. 

In the case of superfluid '^He the SGS is not likely 
to be observable in practice. There are two reasons 
for this. First, the practically achievable weak link di- 
ameters are quite large: generally on the order of the 
zero-temperature coherence length (^o ~ 10 . . . 70 nm), 
and certainly much larger than the Fermi wavelength 
(Af ~ 0.8 nm). Since liquid ■^Hc is naturally free of 
impurities, the quasiparticles simply follow classical bal- 
listics through the aperture. Some non-transparency is 
introduced by scattering at the walls inside a flnite- length 
aperture, but this scattering is diffusive and its principal 
effect is to reduce the net currents Second, also for 
practical reasons, the biases U in weak links of ^He are 
always restricted to the limit U ^ A 12, J^J, while the 
SGS occurs on the scale of A. In fact, in most experi- 
ments U is even much smaller than the quasiparticle re- 
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laxation strength HT due to inelastic scattering, which by 
itself satisfies hT ^ A. The limit of a low-transparency 
point contact between triplet-paired condensates was re- 
cently studied T^l, but, as explained above, such results 
are not likely to be important for interpreting experi- 
ments in superfluid '^He. For intermediate transparen- 
cies, effects similar to those of Ref. may be expected. 

In this paper we consider the limit of a short point con- 
tact with perfect transparency, the so-called "pinhole". 
Furthermore, we concentrate on studying the bias region 
U ^ A, and consider only the B phase of superfluid '^He 
explicitly. However, some of the general results may just 
as well be applied for the A phase, or any other triplet 
or singlet pairing state, and for any value of a constant 
bias U. Even though the SGS in the dc current can- 
not be resolved with our assumptions, there are other 
details introduced by the complicated structure of the 
order parameter in "^He, and its modification due to sur- 
face scattering. The equilibrium limit U = for a •^Hc-B 
pinhole was studied in Ref. ^ in detail, and this paper 
represents an generalization of that calculation to finite 
biases. Parts of our results have already been published 
[lOj l and we review them here in order to obtain a self- 
contained presentation. In addition, we present some new 
analytical results and a more thorough numerical anal- 
ysis of the dc current and supercurrent amplitudes as a 
function of the bias pressure. Some aspects related to the 
so-called "anisotextural" effects are covered in more 
detail elsewhere [l6l| . 

In Sec.ni'WG start with some basic issues of the quasi- 
classical theory, and in Sec. IIIII the pinhole model and 
the general current formulas are introduced. Section 
IIVI presents the analytical results obtained when surface 
pair-breaking is neglected. In this case many limiting 
cases are studied, and we also briefly discuss the con- 
nection of the quasiclassical model to the anisotextural 
effects In. Sec. we present our numerical results 

for the current amplitudes and the sub-gap bound states 
in the presence of the pair- breaking effects, and the com- 
putational methods are briefly explained. A comparison 
of the results to experimental data is provided in Sec. 
IVII and the agreement is found to be good. Section IVlTI 
concludes with some discussion of future directions. Fi- 
nally, details related to the self-consistent computation 
of the order parameters and some mathematical results 
are gathered in the Appendices. 



II. QUASICLASSICAL FRAMEWORK 

Our analysis is based on the nonequilibrium formula- 
tion of quasiclassical theory, which has been throughly 
reviewed in Ref. 0. We start by considering some basic 
points of the formalism here, since it is of essential im- 
portance to the ensuing discussion. The central quantity 
is the Keldysh-space propagator, which has the form 
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(1) 



where g^'^'^ {kjH; e,t) are 4x4 Nambu matrices, and 
"o" denotes the quasiclassical folding product — see 
Appendix^ Here k parametrizes positions on the spher- 
ical Fermi surface of ^He, R is a spatial coordinate, e the 
quasiparticle energy, and t is time. The Nambu matrices 
have the structure 
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gR,A jR,A 
jR,A gR,A 
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where the diagonal components g^'^'^ and off-diagonal 
components f^'^'^ are 2x2 spin matrices, and the con- 
jugation operation " ~ " is defined as g(k, e) = g(— k, — e)*. 
In order to automatically satisfy the normalization con- 
dition in Eq. JQ) , it is convenient to parametrize the prop- 
agator as follows 18, 19, 20]: 
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;^R,A o ^R,A 
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and 

g^ = ^2niN^ 



^R. Q Q ^A-j 



(7^ 



~K 

' O X 



x^ o 7)^ 



X^R Q Q ^A^ 



(4) 



where 



^R,A ^ 



(1 



^R,A Q ^J?..A-) 



(I - ^R,A ^ ^R,Ay 



(5) 

Here the spin matrices 7'^^^(k, R; e, i) are called co- 
herence functions, and they may often be interpreted 
as Andreev-reflection amplitudes. Since they fully 
parametrize g^'^, they completely determine the density 
of quasiparticle states of the system. The spin matrix 
x^ (kjH; e,t), on the other hand, is a distribution func- 
tion describing the occupation of these states. All expec- 
tation values of one-body observables may be computed 
from the Keldysh component g^ , which includes infor- 
mation on the states as well as their occupation. The 
coherence functions satisfy the symmetry t'-'^ = (7"^)^ 
and are related to the spin components of the propaga- 



tor by 7^ — - 
is Hermitian: 
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^ jR The distribution function 



The function x is not the only way to introduce a 
distribution function. A more common definition is given 
by writing 



9 



K 



9" oh 



hog^ 



(6) 



which satisfies the normalization condition for any h. 
However, any physical g^ may be parametrized by choos- 
ing h diagonal 



hi 
-hi 



(7) 
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with the spin matrix hi (k, R; e, t) as the new distribution 
function. Also hi is Hermitian, (/ii)^ = hi, and it is 
connected to by the relation 

oo 

hi = ^ (7-" o7«)"o[a;^-7-«- 0^^07-4] 0(7-4 07^)". (8) 

n=0 

The functions and x^ have the interpretations of dis- 
tribution functions for "particle-like" and "hole-like" ex- 
citations, while hi includes contributions from the coher- 
ent Andreev reflections between the two types. In equi- 
librium hi reduces to the function h{e) = tanh(/3e/2) = 
1 — 2/(e), where /3 = I/ZcbT, T is the temperature, and 
/cb is Boltzmann's constant, and /(e) is the Fermi dis- 
tribution. In comparison, x^ takes the form x^ {e) — 
h{e){l — 7^7'^). Thus hi has a more direct interpre- 
tation as a "quasiparticle" distribution function in, for 
example, the Andreev bound states inside the weak link 
or a vortex core. 

The propagator g satisfies a transport-like equation 
of motion, which depends on self-consistently computed 
self-energies. The latter have a similar Keldysh-space and 
Nambu-space structure as g in Eqs. ^ and P)) fT?! [20l| . 
The equation for g may be rewritten as a Riccati-type 
transport equation for the coherence functions 7^-^, and 
a kinetic equation for x^ |2n| |. In particular, the equation 
for 7^ is 

i^VF • V7'" - 267-^ - 

+ 7«' o A« o 7« + S« o 7« - 7« o E« 

(9) 

where the spin matrices and are the Nambu-space 
diagonal and off-diagonal self energies, respectively, and 
vp = wpk is the Fermi velocity. However, we only need 
these in equilibrium, where the kinetic equation is always 
solved by x^ — h{l — 7^7'^), and the folding products 
in Eq. Q simplify to matrix products. In the mean-field 
approximation E^'^ = E™^ and A-'^'^ = A™^ which 
are independent of e. Most importantly, the off-diagonal 
spin matrix A™^ determines the order parameter of the 
superfluid. In this paper the strong-coupling effects, i.e., 
inelastic quasiparticle-quasiparticle scattering, are only 
taken into account with a simple "normal-state" model 
Y,R,A ^ 5]mf ^ iri(e), where ri(-e) = ri(e). This effec- 
tively adds an imaginary part to energies: e — > e^'^ = 
e ± iri(e). Physically, the imaginary part describes a 
finite quasiparticle lifetime, which is important in the 
parameter ranges of ^He weak link experiments [l^.[T3|. 
Mathematically, it is important for regularizing the diver- 
gences in the MAR process, which occur at low pressure 
biases . In general, the coUisional self-energy gives 
strong-coupling corrections also to A''^'^ (and hence a 
gap-dependent contribution to the lifetime but their 
proper calculation is too complicated for the purposes of 
this paper. 




FIG. 1: Quasiparticles hitting the wall outside of the con- 
striction are scattered, which leads to a suppression of the su- 
perfluid state at distances closer to the wall than (dashed 
trajectories). Only quasiparticles hitting the constriction di- 
rectly are ballistically transmitted and contribute to the cur- 
rent (solid trajectories). 



III. PRESSURE-BIASED PINHOLE 

We now apply the above formalism to describe a small 
constriction of diameter d and area S = 7r(d/2)^ in wall 
between two volumes of superfluid '^He, when there is a 
pressure difference P between the two sides. As a result 
of thermomechanical effects, there may also exist a tem- 
perature difference. Although we initially allow for this 
possibility, we shall always assume the sides to be in good 
thermal contact and thus at equal temperatures. We use 
the so-called "pinhole" model, which is a direct general- 
ization of that used already in Ref. — see Fig. Thus 
we assume that while Ap , it is still much smaller than 
the zero-temperature coherence length ^0 = hv-p 
We also assume that the wall is negligibly thin in compar- 
ison with d, so that scattering inside the aperture need 
not be considered. (A simple way relax the latter as- 
sumption was considered in Ref. |ll|.) The convenience 
of this model is that the calculation of the nonequilib- 
rium current through the constriction need not be done 
self-consistently, since all feedback effects away from the 
aperture may be neglected to lowest order in c?/^o- It is 
enough the know the equilibrium propagators (or coher- 
ence functions) calculated close to the wall on its left 
(1) and right (r) sides in the absence of the constric- 
tion. However, in anisotropically paired superfluid '^He 
the computation of these propagators still requires a one- 
dimensional self-consistent calculation, since t he p resence 
of the surface leads to pair-breaking effects (21| and to 
the existence of surface-bound quasiparticle states below 
the bulk gap. 

In the following we assume that the order parameters 
and the corresponding equilibrium coherence functions 
7-'^* and 7"^* have already been calculated on both sides 
i = 1, r — see Sec. Ivl and Appendix^ for more details. 
Thus the results of this section are still very general and 
applicable to any type of pairing state. 

Let us choose the coordinates such that the z axis is 
perpendicular to the wall and points from 1 (z < 0) to r 
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(z > 0). The current through the pinhole is given by ^3 



de 
27ri 



fc.TrC(e) 5''(k)-.g^(-k) 



(10) 



where (7^(k;e,t) is the diagonal Nambu component of 
inside the pinhole, (■■■)^,^^^ = f^^^^{dfl^^/4'K){- • • ), 
and Tr is a spin-matrix trace. For particle current the 
spin matrix C(e) = 1. For heat current C(e) — e and the 
spin current for spin projection along the axis a — x,y, z 
would be obtained by using the corresponding Pauli ma- 
trix C = (Ja- The unit = \vyN{Q)S = M/nh is the 
normal-state conductance, where A'^(O) is the single-spin 
density of states in the normal state and M is the 
number of conducting transverse modes. 

If we count energies from the chemical potential of the 
side r, the pressure bias P causes a shift in the 1-side 
chemical potential by U = jj} — , where U — [m^/ p)P , 
7713 being the mass of a '^He atom and p the mass den- 
sity of the liquid. The phase difference (f) — ^p^ — 1^9' 
between the 1 and r condensates then varies according to 
(f) = 2U/h. Assuming U to be constant, this is solved by 
(/)(i) — LOjt where the Josephson frequency uj — 2U /h. 
The constant bias U gives 7^' a simple time dependence, 
7-^'(e) 7^\e,t) = e~''^(*^7^'(e), while remains 
time-independent. This allows one to evaluate the fold- 
ing products in [Eq. Ql] analytically — see Appendix 
ini First it is convenient to define 



I{t)^Gn{kl{Kt) 



fe^>0 



(11) 



where /(k, t) is a "channel-resolved" current. Since this 
is periodic with the Josephson period Tj = 211 /ujj, we 
expand 



/(k,i)= ^ /™(k)e- 



-\mu}jt 



(12) 



where /^(k) = /*„(k) = (/(k, i)e™"'^*)T,, and we de- 
fined the time average (• • • )tj = Tj~^ J^' (• • • )dt. Eval- 
uating the folding products one finds, for m > and 

fc. >0 



00 

lm{k)=-J2 I deTrC(e) 

n=0 ' 



n+m,n 



k, e- {2n + m)U, U] 



(13) 



F' 



in,n+m r 



-k,e+(2n + m + l)[/, -[/] 



where 



xi^^{e-Un^\e)][PUe,U)]^ 



U)^Y[ T^'Xe + 2pC/)7^^(e + {2p - l)U). (15) 

The distribution function is x^* — h'^{e)[l~^^^^^^]^ and 
again i = 1, r. The k dependences have been dropped for 
clarity. 

If C(e) is assumed to be energy-independent, then 
Eq. (|13() can be simplified by changing integration vari- 
ables. In the following we also assume the two sides to 
be at equal temperatures, such that /i''''(e) = /i(e) = 
tanh(/3e/2). When the normal-state contribution pro- 
portional to / de[ft.(e + [/) — ft.(e)] — 2U is separated, one 
finds 



/™(k) =TrC<^ 2U5, 



E 

n=0 



[F,';+™'"(k,e,[/)-F, 



de 



(-k, £,-[/)] 



(16) 



Identifying the coherence function 7^ as an Andreev- 
refiection amplitude, Eq. H16I) has a clear interpretation 
as describing the MAR process, with the index n running 
over the number of successive reflections. The present 
results have been derived by assuming U to be constant. 
However, it may be shown that even when U varies in 
time, corrections to the results are small at least if [/ ^ 
hVl. When U < hVi and U < huj^, this should be well 
satisfied. 

The current may also be Fourier-expanded as 



liU, t) =/o(C/) + ^ [liiU) sm{mu;,,t) 



m—1 



(17) 



+ I^{U)cos{mujjt)] 



where the coefficients /o,/^j,/^j are real- valued. They 
are connected to the complex amplitudes of Eq. (|12|) by 



Jo = Gn(fcz/o(k))^^^p, 

^2G„Im(fc,/„(k))^^>„, 
I^,^2GnRe{kJm{k))i^^^,. 

From Eqs. Q and ||SJ) we also note that {g^)Tj 
and thus the dc component is given by 



/o = -G„ / de ( fc, Tr G(e) hi (k) - hi (-k) 



(18) 



~27ri/ii 



(19) 

As seen in Eq. (|16|) . it is convenient to separate the dc 
current as Io{U) = GnU + Iar{U). Here GnU is the 
normal-state part and and Iar{U) is due to MAR only 
At high biases Iar saturates, and gives rise to the 
"excess current" on top of G„J7 (see below). In what fol- 
lows we shall be interested in calculating Iq, ifmlm both 
analytically and numerically for the case of superfluid 
■^He-B. We only concentrate on analyzing the particle (or 
mass) current, where G = 1. 
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IV. RESULTS FOR THE CASE WITH NO GAP 
SUPPRESSION 

A. General results 

For simplicity we shall first neglect the suppressing ef- 
fect of the solid wall on the order parameter. This makes 
the problem formally similar to the s-wave case, and the 
results of this section are rather straightforward general- 
izations of those of Refs. 

i, II II i- The B-phase order 
parameter ij is thus assumed to be of the form A*(k, z) 
= A'(k) • 0-1(72 for j = 1, r, and S™^ = 0. In this case Eq. 
Q is easily solved |0. The gap vectors for momentum 
direction k are given by A'''^(k) = Ae"^' 'd''''(k), where 
d^-'{k) = i?'^'k. Here R^-' = R{h^-\eL) are rotation ma- 
trices, with 9l — arccos(— 1/4) the dipole-locked rotation 
angle 0, and n* the rotation axis on side i = 1, r. 

If, for each k, we choose the spin quantization axis 
parallel to x d"^, the condensates may be divided into 
tt and II parts, which behave much like two independent 
s-wave systems p^ . If we define as the azimuthal 

angles of d''"' in the plane perpendicular to d' x d'', which 



satisfy 



-I- TT, then A and 7 are diagonal 



A^^^ A 



(20) 



(21) 



Here 7^'^ = -A/(e^'^ ii^/A^ - (e«'^)2), and e^-^ = 
e±iri(e), where Fi is present to model inelastic scatter- 
ing. The phase differences of the two condensates over 
the contact are given hy (f)^ ^ = (j) — <^X^i where a — ±1, 
(j) = ip^ — = ujjt, and Xk = arccos(d' • d''). Using these 
definitions, Eq. (|15|l simplifies to [lOl | 



/™(k)=TrC<^2[/,5„,o + 2 







(22) 



xV J detanh(/3e/2)(l - \Y {e)r) 

n=Oq=l p=l+l ^ 



The distribution function hi is proportional to the unit 
matrix in spin space. For fc^ > it is given by 

00 n 

^r(e)-Me) + Eni7''(^-^'^)i' ,,,, 

ri=Oj=0 (23) 

X [h{e -{n+ l)U) - h{e - nU)] 



and for fc^ < by 

00 n 

h<{e) =h{e + f/) + E n It""!^ + U + 

X [h{e + (n + 2)U) - h{e + (n + l)nU)]. 

The dc component of the particle current [Eq. p9|) ] may 
now be written 

/o = -(Gn/2) j de[h>{e) - h<{e)]. (25) 

We note that the result for /n is exactly the same as 
for an s-wave superconductor [a,l3- In particular, it is 
independent of the spin-orbit rotation matrices (iflj . 

B. Limiting cases 

In the small-bias (or adiabatic) limit C/ <C A the vari- 
ation of phase difference (pit) = ujjt is slow. In this case 
one may describe the junction in terms of the occupation 
of the Andrecv bound states |0 

= - Sign(fc, sin((/.^^^/2)) A cos(0t^^/2), (26) 

which are obtained from the poles of in equilibrium. 
The Keldysh function may be now approximated with 
the "quasi-equilibrium" form = hf{g^ — g^) so that 

(27) 



Defining the bound-state occupation probabilities 

Pk,. = {l-/^^k,.W]}/2, (28) 



the particle current may be written 

a.5=±l \ ^ 



(29) 



Neglecting Andreev reflections for |e| > A, we may ap- 
proximate 7-^ w -e-'''('^)-';(')6l(A - |e|) where d{e) = 
arccos(e/A) and C(e) — Fi(e)/vA^ — ? which is strictly 
valid only for |e| <C A. Using these we may approximate 
as 



h>ie)^hie)-eiA-\e\) f de'h'{e') 

J -A 

xexp[-|£de"C(6")] 



(30) 



and a similar expression exists for hf. Then it may be 
shown that the occupation probabilities satisfy the ki- 
netic equation 



(31) 
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where / = {1 — h)/2 is the Ferrai distribution, e = 
e^^[4>{t)], and we defined the relaxation rate T = 
2Ti/h. In a normal-state Fermi-liquid approximation 
r ^ {nk-QTY + , while in the superfluid state some cor- 
rections from the existence of the gap may be expected 
0. The initial conditions for this equation are mostly 
determined by the "thermalization" of the bound states 
when they hit the gap edges at e = ±A 6, J^OJ. Thus if 
at i = to we have ej^ ^.['/'(io)] = ~A (for kz > 0), then the 
occupation is returned to equilibrium: p^^ ^ito) = /(^^)- 
In the limit U <^ hT Eq. may be solved to yield 




hr{e) cosh2(/3e/2) 



(32) 



where e = e^j^. ^.[(/'(t)] and Is{4') is the supercurrent of Ref. 
|2^ Approximating r(e) sa r(0) = Fq, its time average 
may be simplified to 



Io{U) = {A/nro)g{T)GnU. 



(33) 



where the temperature factor is given by g{T) — 
J^-^ tanh(/3Aa;/2)(a;/\/l — x'^)dx. This result is only cor- 
rect to leading order in ?iFo/A, since we have neglected 
the corrections from energies |e| > A. This approxima- 
tion is valid for temperatures not too close to the critical 
temperature Tc. Another exactly soluble limit of Eq. H31(l 
is that of <C A at zero temperature, if we additionally 
assume the Fermi-liquid form r(e) = ce^. In this case 



Io{U) = 2G„A / dxexp 



cA2 
"2[7 



(arcsinx — xyl 



(34) 

which varies as Iq{U) ~ U when U ^ 0, and not lin- 
early as Eq. H33|l at higher temperatures. This is plotted 
in Fig. |2] for experimentally feasible parameters. 

For U ^ hr quasiparticles in the MAR cycle begin to 
reach the gap without being scattered, and Iar begins 
to saturate. In the limit HT <^ U <^ A the current is 
slowly varying, and on the order of 2GnA. Finally, in 
the large-bias limit U ^ A one finds that Iar saturates 
with the asymptotic behavior Q 
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G„tanh(/3C//2). 



(35) 



This is known as an "excess current" , since the total dc 
current is then of the form Iq{U) — G^U + /excj where 
the first term is the normal-state value and the second 
one approaches a constant. 



textures to be fixed. Thus also Xk must be constant 
in time. However, there are situations where the tex- 
tures may also oscillate resonantly with the Josephson 
freauencv|ll|. and the theory may be generalized to take 
such "anisotextural" effects into account. Consider in 
particular the limit U <^ HTq and T ^ Tc, which is real- 
ized in the experiments of Ref. 13, for example. Starting 
from Eqs. ()29|l and (|31|l we find the following expression 



I{t) =7rG„/?A2 / fc^ I cos Xk si 



U 



COS Xk COS (j) + sin Xk sin ( 



1 



(36) 



fc^>0 



where — Lujt and Xk — Xk(^)- The first and third terms 
correspond to the m — 1 terms in Eq. (|17|l . while the sec- 
ond term is the equivalent of Eq. (|33l) . The last term is 
new, and it is only present when Xk is time-dependent. 
Assuming now that Xk(0 oscillates Tj-periodically, then 
we see that all the (^-dependent terms can also have fi- 
nite time averages. Actually, the averages of the last two 
terms in Eq. H36(l exactly cancel each other, but the av- 
erage of the first term gives a dc current in addition to 
Eq. I|33|) . Thus, for T « Tc we may write 



/dc,totai(f/) ~Io{U) + (/f ({/, t) siniLUjt)), 



(37) 



where if has the Tj-periodicity of the angles Xk(i)- The 
first term in Eq. H37|) results from the MAR process, 
where energy is dissipated directly to the quasiparticle 
system. The second term corresponds to dissipation via 
the excitation of collective order-parameter modes, i.e., 
spin waves, which are driven by the oscillating Josephson 
spin currents The true magnitude of the resulting 

dc current depends on the details of the process, which 
is in general geometry-dependent. Therefore, we shall 
not attempt to explore this issue any further here. Note, 
however, that while the dc current Io{U) is independent 
of B}'^ , the second contribution may depend strongly on 
them. In fact, for U hV the coefficient if is generally 
of the form lf{t) = {2m3/h)[aiiizz + a2{ipxx + '4'vv)]{t), 
where tpij = B^^^R^^y The ?7- independent parameters 
ai(T) and a2{T) may be computed numerically even 
when gap suppression effects are accounted for [llj, and 
expect Eq. H37|) to remain valid also in this case. Finally 
note that very similar dc current contributions may arise 
from the "Shapiro" or "Fiske" type effects, where loj co- 
incides with some resonance frequency of the cell instead 
of the order-parameter texture 



quent 

m 



V. GAP SUPPRESSION AND NUMERICS 



C. Anisotextural effects 

Most of the above results have been derived by assum- 
ing the bias U to be time-independent, and the spin-orbit 



A. Calculation of the coherence functions 

For a more realistic calculation, one must take into 
account the suppression of the p-wave order parameter 
Aq = Ao • (Ti(T2 close to solid walls. However, when 
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this effect is included, the calculation of the coherence 
functions may only be done numerically, and the current 
must be computed from Eq. Hlfifl . 

We assume the surface to have translation symmetry 
in its plane, and a rotation symmetry around its normal 
z. If we also neglect external magnetic fields and flow 
parallel to the surface, then the B-phase order parameter 
may be parameterized by the form 



Ao(k,2) = [A||(z)fc„A||(z)fcy,A^(z)fc,], 



(38) 



where z = is at the wall. The real amplitudes A||(z) 
and A^(z) both approach the value A as \z\ oo. Close 
to the wall their behavior must be calculated using one 
of the different models for the surface scattering, which 
are generally ex pres sed as a boundary condition for the 
propagators [2J, 123, l2a, 123 ■ The simplest surface model 
assumes a completely specular scattering of the quasipar- 
ticles from the surface. This already leads to significantly 
improved results for the pinhole currents 11]. However, 
most surfaces are believed to be microscopically rough, so 
that they scatter quasiparticles diffusely. There are sev- 
eral model boundary conditions for such surfaces, which 
all yield practically the same profiles A||(z) and A^(z) 
HI El III El. As in Ref. [Ill we use perhaps the sim- 
plest one, the "randomly oriented mirror" (ROM) model 
|27|. The gap-suppression effect also introduces spon- 
taneous spin currents which flow parallel to the surface 
psj l . This is responsible for the existence of a small spin- 
vector part in the mean-field self energy E™^ = i>'o • cr, 
which depends on the Fermi-liquid parameters with 
odd I |11|. The self-consistency procedure is simplest to 
do by using the imaginary Matsubara energies, as briefly 
discussed in Appendix^ In this part of the calculation, 
the inelastic relaxation rate T is taken to be infinitesimal. 

Once the mean-field self-energies Aq and E™^ have 
been self-consistently computed (see Appendix . the 



coherence functions at real energies 7g 
by integrating the Riccati equation 



R,A 



are calculated 



2e^^i - Ao - Y^'^U 



ymf R, _ ^liyvi 



(39) 



on several trajectories R — uk passing through the pin- 
hole at M = 0. To find the propagator at u = 0, 
only integrations from the bulk toward the wall are 
needed. In Eq. ^ we again introduce the quasiparti- 



cle relaxation rate F = 2Ti/h through e 



i(ftF/2). 



We use the simple normal-state Fermi-liquid form F = 
a[(7rfcBr)^ + e2]/(ro7r2A:|), where tq = 1.14 us mK^ is 
obtained from viscosity measurements [H E^ , and a is 
a free parameter of order unity. The bulk solution of Eq. 
(|39|l is easily found El], and this is used as an initial 
condition. For a junction with mirror-symmetry with 
respect to z = 0, the 1 and r solutions at the junction 
satisfy 7^''(-k, z = 0) = -[7^'(k, z = 0)]^, and thus the 
integration need only be done on one side. The different 
spin-orbit rotations of the B phase on the two sides of 




0.02 



^ B c ' 



FIG. 2: The dc currents la for temperatures T/Tc = 
0.4, 0.5, 0.6, 0.7, 0.8, 0.9 in order of decreasing amplitude. The 
solid lines include the gap suppression effects, whereas in the 
dashed lines it is neglected. The dash-dotted line is the zero- 
temperature result of Eq. 1341 , while the straight dashed line 
corresponds to Eq. 1331 at T/Tc = 0.4. The results are simi- 
lar for both parallel and antiparallel n^'^^'s. Here Fi — and 
a = 1.6. 



the junction may be taken into account with the trans- 
formations 

7«^(k,z,e) = C/^f (k,z,e)rF, 

A^k, z) = WAlik, z)[uY, (40) 
I]'"f'*(k,z) = U'J:f'\-k, z)[W]\ 

where W — exp(— i6'Ln' • <t/2), and i = \,r. Finally, the 
functions 7^'(z = 0) are inserted into Eq. Hl()|) to obtain 
the current amplitudes. 

We note that the introduction of F in Eq. H39|l is not 
self-consistent, and stress that our normal-state model 
neglects all strong-coupling modifications on A^'"^ [l7l| . 
This procedure should therefore be regarded only as a 
rough model for the description of the quasiparticle- 
quasiparticle scattering. However, the calculation still 
provides at least a semi-quantitative model for studying 
the simultaneous effects of gap suppression and inelastic 
processes on the pinhole currents. 



B. Current amplitudes 

Figures f2l6l show our results for the lowest-order cur- 
rent amplitudes Iq, 1^2 and ^-s a function of the bias 
U. The amplitudes and for m > 3 may be set 
to zero since they are negligibly small for all practical 
purposes. Figures El and ^ are for parallel n's on the two 
sides of the contact (n' = n'' = z) and Figs. [S] and |H| 
are for the antiparallel case (— n' = n'' = z). The Iq 
amplitude shown in Fig. |21 is the same for both cases. 
Actually, as our numerical calculation shows, the result 
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FIG. 3: The amplitudes if and if for T/Tc = 0.4,0.6,0.8 in 
order of decreasing amplitude in the case of parallel n'''^'s. The 
solid and dashed lines are for if with and without gap sup- 
pression, respectively, and the dash-dotted and dotted lines 
are the corresponding values for /p. Here _Ff — and a = 1.6. 
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FIG. 4: Same as Fig. El but for if and I^ 



for /q is practically independent . For the case where 
gap suppression is neglected, this independence was pre- 
viously shown to be exact — see Eqs. and H25|) . 

In the calculations, we have always assumed — 
for Z > 3, while the effect of the remaining with values 
between —1 ... on the pinhole currents is at most a few 
percent of their total amplitude — in the figures of 
this paper = also. The quasiparticle relaxation 
parameter is chosen as a = 1.6, since this is the value 
which gives the best fit to the relevant experiments — see 
Sec. lVIl Although the behavior of the current amplitudes 
for U < hT depends strongly on the choice of a, the 
asymptotic behavior for U ^ HT does not. 

We have only calculated the currents numerically down 
to T/Tc = 0.4, although Fig.|21shows the additional zero- 
temperature result of Eq. (|34|l . This is because the re- 
quired number of terms in the MAR sum [Eq. (|15|l ] is 
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FIG. 5: Same as Fig.Elbut for antiparallel ^''^^'s 
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FIG. 6: Same as Fig. gl but for if and I^ 



proportional to A/Fq, which is on the order of thou- 
sands already at T/Tc — 0.4. We have taken into ac- 
count all terms up to n = 4A/?iFo. Actually, this is 
only necessary for U <^ hV where the number of suc- 
cessive Andreev reflections is limited by relaxation. In 
this regime it would also be possible to use approximate 
schemes instead of the full expression (|l(j|l . On the other 
hand, for U > HT a maximum of ^ A/U terms should 
be enough. The energy cutoff for the coherence functions 
was typically chosen at around lOfceTc, and the density 
of energy discretization points was highest close to the 
gap edges, where the highest accuracy is needed. For en- 
ergies in between these points, linear interpolation was 
used in computing the current amplitudes with Eq. 
The number of Gaussian polar angles used in angular 
integrals was usually eight. 

In the limit [7^0 the cosine amplitudes all van- 
ish (hnearly in U/HTq) as they should, since the equi- 
librium current Is {(f) must satisfy the time- reversal sym- 
metry Is{—(j>) = —Is{(j>). Upon inserting the ampli- 
tudes to Eq. (|17|l in this limit, the current-phase relations 
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a=+1 




FIG. 7: Local density of states inside the pinhole for kz — 
0.93 and — n' — A'^ = z when gap suppression at a diffusive 
wall is included. Here a large, arbitrary value a ~ 75 was 
chosen for purposes of illustration, in order to make the width 
~ hr of the bound-state peaks better observable. Among the 
peaks, those corresponding to the a = +1 {a — —1) spin 
branch are the ones apparently shifted toward smaller (larger) 
0. Note that the slope of the bound states is not so steep as 
for Eq. 12611 . Also, for given phase difference and spin band, 
more than one bound state can exist simultaneously. 



Is{(j)) of Ref. are quite accurately reproduced, when 
we make the replacement ujjt — > (jj. Note, in particular, 
that in the antiparallel case with no gap suppression, the 
second-order amplitudes are of equal magnitude with the 
first-order ones. Thus Is{4>) exhibits a strong "tt state" 
[2^ . When the gap-suppression effect is taken into ac- 
count, the m = 2 amplitudes tend to be very strongly 
suppressed, and in this case the "tt state" only presents 
itself at very low temperatures Thus as a first ap- 
proximation, j| may often be neglected in comparison 
with if. However, apart from the limit U HT oi very 
low biases, the cosine amplitude /f as well as the dc 
component Iq are equally large as the sine amplitude if . 
Therefore the models which are based on if alone — like 
that of Ref. ^3 — are only valid for U HT and T close 
toTc. 



density) 

iV(k,R;e) = -^^ImTr5^(k,R;e). (41) 

TT 

These peaks correspond to poles of g^, which depend on 
the phase difference (p — ip'-' — ip^ and the rotation matrices 
through V'y = B}^^R^^y The width of the peaks is on the 
order of the relaxation rate hT. Although the spectral 
weight of the bound states is spread a round the pinhole 
up to distances of order hv-p / a/A^ — e^, we only calculate 
Eq. igTl) inside the hole, at R = 0. 

Figurc[71illustrates the results for one choice of ^jjij and 
k with = . . . 27r when the gap suppression effect of a 
diffusive surface is included. Since ^ 5ij, the bound- 
state peaks in the densities show clearly a spin-splitting 
between the "ct = ±1" condensates, as in the simple an- 
alytic result (|26|l . Since the splitting angle depends on k, 
an average of A^(k, e) over the Fermi surface in fact leads 
to the a formation of a wide hand of bound state ener- 
gies. The bound states for — k are obtained by using the 
time-reversal symmetry k — s- — k, <j) — > —0, a —a, and 
the "particle-hole" symmetry k — k, e — e, which 
follows from the symmetries of the equations of motion 
and the geometry |25l |. 

Compared to Eq. it is seen that the gap sup- 

pression modifies the bound states such that they are 
always at energies |e| < A. In the bulk A^(k, e) « 
2iV(0)6'(|e| - A)|e|/Ve2 - A^ has divergences at |e| = A, 
but in the middle of the junction most of the spectral 
weight is now in the bound states even at = 0. Thus 
there are several branches of bound states coexisting si- 
multaneously for given (p. Accordingly, Eq. (|29|l for the 
current should be modified by replacing with 
e ^^^(0), and by adding a sum over the branch index 
q. Since the energies e^j^^(0) for a given branch in the 
range (—A, A) are now mapped to phase differences in 
the range (—00, 00), their slopes are not so steep. There- 
fore, one would expect that the dc currents are generally 
smaller when the gap suppression effect is taken into ac- 
count. As seen in Fig. |21 this is usually the case. 



VI. COMPARISON TO EXPERIMENT 



C. Bound states 

In Sec. II VI we saw that in the adiabatic regime U ^ A, 
the currents in the point contact may simply be described 
in terms of the equilibrium bound states. It should be 
possible to generalize this description to the case where 
gap-suppression is included. To study this, we have cal- 
culated numerically the density of quasiparticle states at 
the pinhole. The bound states are given by the sub-gap 
peaks in the k-resolved local density of states (or spectral 



In this section we present a brief comparison of the 
above theory to available experimental data on the 
current-pressure characteristics in ■^He-B weak links 0, 
H^ . There are a couple of basic things to note about 
the experiments. First, as already mentioned, it is dif- 
ficult to manufacture apertures which would satisfy the 
requirements of a pinhole very well The apertures 
with d « 100 nm in a 50 nm membrane used in Ref. [3 
are rather close, at least compared with the 0.25 /.tm wide 
slits in a 0.1 fim membrane of Ref. [l^. Second, in order 
to ensure leak-proofness, the pressure biases are limited 
to very low values where i7 <C A. The experiments of 
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FIG. 8: Comparison of the data (+ signs) from R ef . to the 
pinhole theory (sohd hnes) for a diffusive surface and a « 0.27 
at indicated temperatures. 



Ref. |l3 are even restricted io U <C fi-F ^ A. There- 
fore, the "excess current" hmit of Eq. (|35|l . for example, 
seems not practically achievable in superfluid '^He. Nev- 
ertheless, the data of Refs. 111113 are enough to make a 
comparison between the most important features of the 
theoretical and experimental current-pressure character- 
istics. 

Figure IHl shows a comparison between the data of Ref. 
Il^ and the numerical pinhole calculation of Io{U). A 
diffusive surface is assumed, and we use the parameter 
a K, 0.27. Again we note that although the slope at 
U — Q depends strongly on a, the asymptotic behavior 
for U ^ hV does not. Indeed, at J7 > ftF the agreement 
is rather good for any value of a on the order of unity, 
although a perfect fit for all temperatures simultaneously 
cannot be achieved. Note also that the experimental cur- 
rents are not even approaching zero in the limit C/ — > 0. 
However, a better fit can hardly be expected, since the 
apertures used in these experiments were far from good 
pinholes. In fact, in sufficiently large apertures one would 
expect a transition into a regime where the dissipation is 
best described in terms of phase slips by vortices rather 
than MAR. Studying this transition would be interesting, 
but computationally very demanding. In any case, the 
overall form of Iq{U) is very similar to the experimental 
results, and the order-of-magnitude agreement on both 
axes is surprisingly good. This gives strong support to 
the expectation that the dominating dissipation mech- 
anism in small apertures of superfluid ■^He is the MAR 
process, analogously to superconductor point contacts. 

The apertures used in the experiments reported in 
Ref. llSlare somewhat better approximations to pinholes. 
Therefore we have chosen to use these data to estimate 
the value of a for the numerical calculations of the pre- 
vious Section. The experiments were carried out in the 
limit U ^ hr, and according to Eq. H33|) . the current 




20 

pressure (mPa) 



FIG. 9: Comparison of the L state data {+ signs) from Ref. 

to the pinhole theory (solid lines) for a diffusive surface 
and a ~ 1.6 at indicated temperatures. 



should be linear in the bias. For the "L state" this is 
rather well satisfied, and a fit to the L-state data gives 
a « 1.6, when gap sup pression at a diffusive surface is 
taken into account . The fit is shown in Fig. |51 The 
increasing low-temperature deviations in this fit may be 
partly due to the insufhciency of our normal-state model 
for F. The currents for the "H state" reported in Ref. 

are larger and more nonlinear than in the L state. 
We suspect that part of this nonlinearity may be due to 
the additional dissipation effects described by Eq. l|T7| 
|l(i| . Another possibility is that the experimental aper- 
tures already deviate so strongly from pinholes, that the 
dependence of the bound states on textures is not suffi- 
ciently described by simple phase shifts. 



VII. CONCLUSIONS AND DISCUSSION 

In conclusion, we have presented an analysis of 
pressure-biased weak links between two volumes of su- 
perfluid ■^He-B by using the pinhole model of a short, bal- 
listic point contact. We showed how the s-wave results of 
Refs.0,|3may easily be reproduced and generalized to the 
p-wave case by parametrizing the Green functions with 
the so-called coherence functions. In the case where the 
gap suppression at surfaces is neglected, we calculated 
the current amplitudes analytically in several limits. In 
the general case, the order parameters and the pinhole 
currents were calculated numerically. Comparison to ex- 
periments gives strong support to the existence of the 
MAR effect in pressure-biased weak links of superfluid 
■^He P^]. We also predict the existence of additional dc 
current contributions, which result from the excitation of 
collective order-parameter modes. These "anisotextural" 
effects are discussed elsewhere in more detail [T6| . 
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In order to improve upon the results of the present 
paper, one should take into account the strong-coupling 
effects in a more detailed way than with our "normal- 
state" model for the quasiparticle relaxation rate. One 
should also consider apertures of finite size, at least by 
computing the bound-state spectra in equilibrium to see 
if textures have some significant effect on them in this 
case. A dynamical calculation for the finite-size aper- 
ture should also be performed, but this already seems 
to approach the limits of practical feasibility. A fully 
self-consistent calculation of the anisotextural effects in 
a pressure-biased pinhole array also appears to be very 
difficult. However, until such improvements are made, 
a parameter-free comparison with the experimental data 
cannot be expected. 

With somewhat less effort, the pinhole heat conduc- 
tivity 30] or spin currents could be studied by starting 
from Eq. (|13|l . The current-noise properties Q of a '^He 
pinhole could also be of some interest, for exarnple in the 
design of accurate superfluid "^Hc gyrometers jSlj. 
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APPENDIX A: EQUILIBRIUM EQUATIONS 

Close to the planar wall at z = the mean-field self en- 
ergies and the coherence functions must be iterated self- 
consistently. Since we are concerned with an equilibrium 
system, it is easiest to do this by using the Matsubara 
technique, where one makes an analytical continuation 
from the real energies e^'^ = e ± iO"*" to the imaginary 
Matsubara energies ibi|em|, where em{T) = 7rfcBT(2m+l) 
and m = 0,±1,±2,... [l3|. The coherence functions 
are then obtained from Eq. I|39|) . where —> ie„i and 
7(f (^) loi^m)- (Here we take Fi as infinitesimal.) 
The self-energies are of the form Aq = Aq • <Ti(T2 and 
E™ ^ = Vq ■ (T. If we use the symmetry i>o = — i/q 
[nl [13, the solution of Eq. ^ may be parametrized 
by 7o = To ■ "'io'i/ 1 a-nd the equation becomes 



where the p-wave pairing interaction y(k • k) — SVik • k' 
was assumed, and 



i/iVF • V7o = - 2ie„7o - Aq - 2(7o • Ao)7o 
+ (7o • 7o)^o + 2ii^o X 7o. 



(Al) 



Due to translational symmetry along the wall, the solu- 
tion and the self-energies only depend on the coordinate 
z. The self-consistency equations are given by [17| 



Ao(k,z) = 3fcBTFi J2 ((k-k')fo(k',2,em) 



z^o(k,z) = 3fcBT 

|e™|<e 



Aa(k-k')go(k',z,e„) 



(A3) 



In these (•••)k = /(dr!k/47r)(- • • ), go = Tr2(cr5o)/2, 
fo = — Tr2((T(7oi(T2)/2 where the upper diagonal and off- 
diagonal Nambu components of the Matsubara propa- 
gator 5(k, z,e„) are given by go = -7ri(l + 7o7o)(l - 
7o7o)"^ and /q = -27ri7o(l - 7o7o)"^- The conjuga- 
tion symbol " ~ " now means 7o(k, e„i) — 7o(— k, £„)* 
~ [70 (k, — em)]t. The coupling constant Vi may be elim- 
inated in the usual way by noting that at T ^ the 
gap vector Ao is small and fo ~ 7rAo/|em|- Thus Ao 
may be canceled and one finds 



yf 1 « nkBT, 



E I 

'-m\<ec 



(A4) 



|£™|<ec 



(A2) 



where = em{T). Inserting this into Eq. 1A2|1 . the cut- 
off Cc may be taken to infinity — see Eq. (12) in Ref. ITlI 
The function Aa{x) = ^2^0 ^tl'^ + ^iVi^l + ^T'Pii^)^ 
where Pi{x) are the Legendre polynomials. From sym- 
metries of the planar wall geometry it follows that if k is 
in the xz plane, so are fg, Ao(k) — A||fc2;X + Aj_kzZ and 

7g, while go and i-'o(k) = i'oy(k)z are perpendicular to 
it. Also due to symmetries, all F° with even I drop out 
of the theory, and for ^ > 3 we assume them to be zero 

El 

In the bulk Ao is constant and I'o = 0, so that Eq. 
HA1|) is easily solved, and this solution is used as an ini- 
tial condition on trajectories starting from the bulk. On 
the wall side, one needs a boundary condition to compute 
the "outgoing" coherence functions from the "incoming" 
functions. We use the ROM boundary condition, which 
is explained in Ref.l2^ This is easiest to express in terms 
of the components of the propagator g(k, z,em), rather 
than directly with the coherence functions 7o(k, z,em). 
The relation 70 = — (Itt — 50)^^/0 is then useful for com- 
puting the initial condition for 70 on the outgoing trajec- 
tories 32] . A specularly reflecting surface is much sim- 
pler to implement, since it only leads to the the continu- 
ity condition 7o(k) — 7o(k), where k = k — 2(k • z)z is 
the specularly reflected direction. As shown in Ref. 0, 
the pinhole current amplitudes for specular and diffusive 
surfaces usually have only minor (at most a few percent) 
differences. Thus, in practice the implementation of a 
diffusive surface may not be worth the trouble. 



APPENDIX B: FOLDING PRODUCTS 

The quasiclassical folding product between objects of 
the type A{e,t) is defined in Ref. 0. This result is 
generalized to multiple products of n objects Aj{e,t), 
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j = 1, 2, . . . , n as follows 
(ii o...oi„)(e,i) 

= n + 1(- ■ • + ^^-1 - ^^+1 -■••)' ^^i' ^^^^ 

assuming all the Fourier-transformations leading to the 
"mixed" representations Aj{e,t) exist. Here dj refers 



to a derivative with respect to the jth time variable tj. 
When the time dependence of Aj(e,t) is harmonic [i.e., 
cx exp(iei/?i)], Eq. (|B1|I yields the folding product ana- 
lytically. In equilibrium all time derivatives vanish, and 
the folding product becomes a simple matrix product. 
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